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Abstract. A brief review to string and parton percolation is presented. After a short introduction, the 
main consequences of percolation of color sources on the following observables in A-A collisions: J/tp 
suppression, saturation of the multiplicity, dependence on the centrality of the transverse momentum 
fluctuations, Cronin effect and transverse momentum distributions, strength of the two and three body 
Bose-Einstein correlations and forward-backward multiplicity correlations, are presented. The behaviour 
of all of them can be naturally explained by the clustering of color sources and the dependence of the 
fluctuations of the number of these clusters on the density. 

PACS. 2 5.75.-q, 12.38.Mh, 24.85.+p 
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1 Introduction 

What conditions are necessary in the pre-equilibrium stage 
to achieve deconfinement and perhaps subsequent quark- 
gluon plasma formation? This question on the occurence 
of color deconfinement in nuclear collisions without assum- 
ing prior equilibration has been addressed on the basic of 
two closely related concepts, string or parton percolation 
0-0 and parton saturation 0-0 -[I] .In this paper we 
will study the first subject. 

Consider a flat two dimensional surface S (the trans- 
verse nuclear area), on which N small disc of radius ro 
(the transverse partonic or string size) are randomly dis- 
tributed, allowing overlapping. With increasing density 
n = N/ttR 2 (we take here S — ttR 2 ), clusters of increasing 
size appear. The crucial feature is that this cluster forma- 
tion shows critical behaviour: in the limit N — > oo and 
R — > oo with n finite, the cluster size diverges at a certain 
critical density. The percolation threshold is given by 

and its value 1.13 is determined by numerical studies. For 
finite N and R, percolation sets in when the largest cluster 
spans the entire surface from the center to the edge. Be- 
cause of overlap, a considerable fraction of the surface is 
still empty at the percolation point in fact, at the thresh- 
old, only 1 — exp(—rj c ) ~ 2/3 of the surface is covered by 
discs. 

In high energy nuclear collision, the strings or partons 
are originated from the nucleons within the colliding nu- 
clei, therefore their distribution on the transverse area of 
the collision is highly non uniform with more nucleons 
and hence more strings or partons in the center than in 
the edge. In this case the value of r] c becomes higher 



2 Local parton percolation and J/ip 
suppression 

Hard probes, such as quarkonia, probe the medium locally, 
and thus test only if it has reached the percolation point 
and the resulting geometric deconfinement at their loca- 
tion. It is thus necessary to define a more local percolation 
criterium 0. 

As we mentioned before, at the percolating critical 
density, 1 /3 of the surface remains empty. Hence disc den- 
sity in the percolating cluster must be greater than | -2% . 

In fact numerical studies show that percolation sets in 
when the density of partons in the largest cluster reaches 
the critical value 1.72/7r7g, slightly larger than This 
result provides the required local test: if the parton den- 
sity at a certain point in the transverse nuclear collision 
plane has reached this level, the medium there belongs 
to a percolating cluster and hence to a deconfined parton 
condensate. In Fig. 1 the percolation probability and its 
derivative as a function of 77 are shown. 

Let us apply the above idea to J/ tfj suppression in A-A 
collision |Hj. We denote by n s (A) the density of nucleons 
in the transverse plane and by dN q (x, Q 2 )/dy the parton 
distribution functions. (At central rapidity y — 0, we have 
x = kr / y/s , where kr denotes the transverse momentum 
of the partons and thus kr — Q). The local parton perco- 
lation condition is 



n s {A) 



( dN q {x,< 
V dy 



x=Qc/\/s 



1.72 

vol 



(2) 



For a given A-A collision at a fixed centrality and en- 
ergy, the relation (0) determines Qc- 
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1. Percolation probability and its derivative as a function r* 



For P6-P6 collision at a/s = 17.4 GeV, Q c ~ 0.7 GeV. 
The scales of the charmonium states %c and tp' , as de- 
termined by the inverse of their radii calculated in po- 
tential theory, are around 0.6 GeV and 0.5 GeV respec- 
tively therefore the parton condensate can thus resolve 
these states and all x an d tp' states formed inside the per- 
colating cluster disappear. The location is determined by 
the collision density. The first onset of J ftp suppression 
in Pb — Pb collision at SPS should occur at N part ~ 125, 
where the J ftp's due to feed-down from Xc and tp' states in 
the percolating cluster are eliminated. Directly produced 
J ftp's survive because of their smaller radii (leading to a 
scales of 0.9-1.0 GeV) and its dissociation requires more 
central collisions, which lead to a better resolution, i.e. to 
an increase of Q c . For Q c = 1.0 GeV we need N part ~ 320. 
The resolution scale of the direct Jf tp cannot be reached in 
S-U collisions, therefore only one stop pattern suppression 
is obtained for this case. 

For Au-Au collision at RHIC, the increase parton den- 
sity shift the onset of percolation to a higher resolution 
scale, so that from the threshold on, all charmonium states 
including J ftp's are supressed, starting at N part — 90, i.e. 
a single step suppression pattern occurs. 

For the case In-In collisions at SPS energies, the thres- 
hold for directly produced J/ip's is not reached even for 
the most central collisions, and again a single step sup- 
pression patterm is expected. 

A detailed discussion and comparison with experimen- 
tal data can be found in references and EH- 



qq pairs in this field. With growing energy and/or atomic 
number of colliding particles, the number of strings grows 
and they start to overlap, forming clusters. At a critical 
density a macroscopic cluster appears that marks the per- 
colation phase transition. 

The percolation theory governs the geometrical pat- 
tern of the string clustering. Its observable implications, 
however, require introduction of some dynamics to de- 
scribe string interaction, i.e, the behaviour of a cluster 
formed by several overlapping strings. 

It is assumed that a cluster behaves as a single string 

with a higher color field Q n corresponding to the vectorial 

sum of the color charge of each individual Qi string. The 
resulting color field covers the area S n of the cluster. As 

Q n = Y^l Qu an d the individual string colors may be 
oriented in an arbitrary manner respective to one another, 



-> 2 



-» 2 

n Qi 



the average QuQij is zero, and Q 

Knowing the charge color Q n , one can compute the 
particle spectra produced by a single cluster of such color 
charge and area S n using the Schwinger formula. For the 
multiplicity fi n and average p T of particles, < p\ >„, 
produced by a cluster of n strings one finds |11|-|12| 



n S\ 2 ,o\ 

— <p 2 T >i (3) 



where \i\ and < p\ >i are the mean multiplicity and mean 
p\ of particles produced by a single string with a trans- 
verse area Si = nr^. For strings just touching each other 
S n — nSi and hence /j, n = n/xi; < p T > n =< p\ >i as ex- 
pected (simple fragmentation of n independent strings). 
In the opposite case of maximum overlapping, S n = Si 
and therefore fx n = \/n^i,< v\ >n= y/n < Vr >ii so 
that the multiplicity results maximally supressed. Notice 
that a certain conservation rule holds 



Mrs 

n 



< Pt Ml < PT >1 



and also the scaling low 

< Pt >n fUnSn =< Pt >i / PiSi 
In the limit of high density 



(4) 



(5) 



3 String percolation 

Multiparticle production is currently described in terms 
of color strings stretched between partons of the projectile 
and target, which decay into new strings through q—q pro- 
duction and subsequently hadronize to produce observed 
hadrons. Color strings may be viewed as small discs in 
the transverse space, tttq, tq — 0.2 — 0.25 fm, filled with 
the color field created by the colliding partons. Particles 
are produced by the Schwinger mechanisms jlOj emitting 



< nSi/S n >= 



1 — exp(—rj) F 2 (jj) 



,r) = Nsirrl/TrR' 2 



Thus 



H = N s F(ri)in ; < p T >=< P T >i / F (v) 



(6) 



(7) 



The universal scaling law (J3J) is valid for all projectiles 
and targets, different energies and centralities, being in 
reasonable agreement with experimental data |13j . A sim- 
ilar scaling is found in the color glass condensate approach 
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Notice that N$ ~ N co u ~ N^ 3 at central rapidity (at 
the fragmentation region N$ ~ N part ~ iV^). As -F^) ~ 

Therefore, the multiplicity per participant does not de- 
pend on the number of participants, there is saturation. 
Numerical studies show a good agreement with SPS and 
RHIC data [12]. The prediction for central Pb - Pb colli- 
sions (N A = 400) at LHC Qs = 5.5 GeV) is fi/0.5N part = 
8.6 and the total charged multiplicity per unit rapidity (at 
central rapidity) is 1800. This numbers are very similar to 
the ones obtained from Hera data, assuming scaling, using 
PQCD and the BK equation [32|-[n|. 




<p x >! < Pt > 3 < PT > n 



No fluctuations No fluctuations 

Fig. 2. Cluster formation for low, intermediate and high den- 
sity respectively. 
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4 Transverse momentum fluctuations 



The behaviour of the transverse momentum fluctuations 
can be understood as follows: At low densities most of the 
particles are produced by individual strings with the same 
< Pt >i, so the fluctuations must be small. Similarly, at 
large density, above the percolation critical point, there 
is essentially one cluster formed by almost all the strings 
created in the collision and therefore fluctuations are not 
expected either. Indeed, the fluctuations are expected to 
be maximal when the number of different clusters becomes 
larger, just below the percolation critical density (see Fig. 
2). In this case in addition to the normal fluctuations 
around the mean transverse momentum of a single string, 
there are more fluctuations due to the different average 
transverse momentum of each cluster. 

Experimentally, it has been measured the quantity 
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300 350 400 
Number of participants 



Fig. 3. F PT (0/0) vs the number of participants. Experimental 
data from PHENIX at ^fs = 200 GeV. Solid line our results. 



5 Universal transverse momentum 
distributions 



,-, _ ^Data — ^rand Oil i 

_ V< P| > - <PT > 2 

t pt = U> — 

^random < PT > 

( § ) 

where Lo ran dom denotes the corresponding normalized fluc- 
tuations in the case of statistically independent particle 
emission. 

In Fig. 3 our result |18) compared with the experimen- 
tal data is shown. A reasonable agreement is obtained. 
There is an alternative explanation based on the occur- 
rence at RHIC of minijets which will enhance the px fluc- 
tuations. At high centrality, it is well established, the sup- 
pression of high pt particles at RHIC what can explain 
the suppression of fluctuations seen at lower centralities. 
According to this picture, at SPS where the production 
of minijets is negligible, this behaviour is not expected, 
contrary to the expectations of percolation of strings. 

Instead of pr fluctuations, the NA49 Collaboration 
has measured (< n 2 _ > — < n_ > 2 )/ < n_ >) at SPS 
for Pb-Pb collisions, as a function of the centrality of the 
collision, showing a maximum at low centrality and being 
1 at high centrality. This behaviour has nothing to do with 
minijets. On the contrary, it is explained naturally in our 
approach |277] . 



In order to know the transverse momentum distributions 
one needs the fragmentation function f(x,pr) for each 
cluster, and the mean squared transverse momentum dis- 
tribution of the clusters, W(x), which is related to the 
cluster size distribution through Eq. @- For f(x,px) we 
assume the Schwinger formula, f(x,pr) = exp{—p T x), 
used also for the fragmentation of a Lund string PJJj at 
first approximation x is related to the string tension or 
equivalently to the inverse of the mean transverse mo- 
mentum squared. For the weight function W(x) we choose 
the gamma distribution 

W(x) = jl-ijx^expi-jx). (9) 

The reason of this choice is the following: In periph- 
eral heavy ion collisions there is almost not overlapping 
between the formed strings and therefore the cluster size 
distribution is peaked around low values. Most of the clus- 
ters are made of one single string. As the centrality in- 
creases the number of strings grows, so there are more 
and more overlapping among the strings and the cluster 
size distribution is strongly modified, according to Fig. 
4 where three cluster distributions corresponding to three 
different centralities of the collision are shown. Each curve 
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and the normalized pt distribution is 




Fig. 4. Schematic representation of the number of clusters as 
a function of the number of strings of each cluster at three 
different centralities (the solid line corresponds to the most 
peripheral one and the pointed line to the most central one). 



in Fig. 4 can be reproduced by gamma distributions with 
different k values. 

Moreover, the increase of centrality can be seen as a 
transformation of the cluster size distribution of the type 



P(x) 



xP(x) 

<x> 



x k P(x) 



(10) 



This kind of transformation were studied long time 
ago by Jona-Lasinio in connection to the renormalization 
group in probabilistic theory [22]. Actually an increase of 
the centralities is equivalent to a transformation which 
changes cells (single strings) by blocks (clusters) and the 
corresponding variables /ii and < p\ >i of the cells by \i n 
and < p\ > n - These transformations of the type of the 
chain (10) have been used also to study the probability 
associated to events which satisfy some requirements |23) . 

The 7 and k parameters of the gamma distribution are 
related to the mean x and dispersion of the distribution 
through 



< x >- 
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< x 2 > - < x > 2 

< X > 2 



(11) 



We use Eq (7) to take into account the effect of over- 
lapping of strings, and hence f(x, tut) — exp(—p T xF(r])). 
Therefore we obtain 



1 



fc<pij, >1 



r(k) 



dx exp{—p T xF(rj)) 



{^x) k 1 exp(-jx) (12) 



f(pr,y) 



dN (k - l)F(r)) 



dy k < p\ >j 



1 



F{r,)p 2 T 
k<pL>i 



(13) 



The equation (|12fl can be seen as a superposition of 
chaotic color sources (clusters) where 1/k fixes the trans- 
verse momentum fluctuations. At small density rj << 1, 
the strings are isolated and there are not fluctuations, 
k — > oo. When the density increases, there will be some 
overlapping of strings forming clusters, the fluctuations in- 
crease and k decreases. The minimum of k will be reached 
when the fluctuations in the number of strings per clus- 
ter reach its maximum. Above this point, increasing 77, 
these fluctuations decrease and k increases. In the limit, 
when only one cluster of all strings is formed, there are 
not fluctuations and again k — > 00. 

The obtained power-like behaviour (p T )~ k , with an 
exponent k related to some intrinsic fluctuations, is com- 
mon to many apparently different systems, as sociologi- 
cal, biological or informatic ones. Distributions like the 
citations of scientific works, or other complex networks 
[2S]|2ni where the probability P(m) of having a given 
node with m links is described by the free scale power law 
P(m) ~ (m)~ k with k related to the fluctuations in the 
number of links obey the same behaviour. Also, it has been 
shown |27|-|28| that maximization of the non extensive in- 
formation Tsallis entropy leads to the same distribution 
(12). 

The universal behaviour indicates the importance of 
the common features present in those phenomena, namely, 
the cluster structure and the fluctuations in the number 
of objets per cluster. 

^From l|13|l one can calculate 



din f 
dlnpx 



-2F( V ) 



Pt > 



1 



EM p1 



k 



<p^>, 



< P T >lt ' 



(14) 



where "i" refers to the different particle especies. 
As p\ — > this reduces to —1F{j])p T j <p\ 
This behaviour has been confirmed by the PHOBOS 

Collaboration. As < p\ >\ p > <p\>\k > <p\>\k the 

absolute value is larger for pions than for kaons and than 

for protons. 

Now, let us discuss the interplay between low and high 
Pt- One defines the ration Rcp(pt) between central and 
peripheral collisions as 



Rcp(pt) = 



f(p T ,y = 0)/N> 



coll 



f( PT ,y = 0)/N coU 



(15) 



where the distribution in the numerator corresponds to 
higher densities 77' > 77. In the pt — > limit, taking into 
< 1 and that F(r}') < F(rj), we 



account that | 
obtain 



< 



fe-i 



Rcp(0) 



(16) 
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approximately independent of k and k' . As rf /i] increases, 
the ratio Rcp(0) decreases, in agreement with experimen- 
tal data. 

As pt increases we have 



Rcp{pt) 



1_+ F{r))"p\l < p\ >ii 
1 



F(v')Pt/ < Pt >u 

and Rcp(pt) increases with px (again, F(r)) > F(rj'j) 
At large px, 



(17) 



Rcp(pt) 



Fjrj) k' 
F(i) k 



2 \k-k' 



(Pt) 



(18) 



which means that if we are in the low density regime k > k' 
and Rcp{pt) > 1, and we reproduce the Cronin effect. 
As we increase the energy, the density increases and on 
reaches the high density regime where k' < k and sup- 
pression of occurs. The Cronin effect disappears at high 
energies and/or densities. The critical density at which the 
Cronin effect disappears is the same at which the trans- 
verse momentum fluctuations presents a maximum. 

Rcp(pt) for two different particles, for instance p and 
7r, becomes, at intermediate pt, 



R P CP {Pt) 

Rcp(pt) 



2 \ k'- 

< P T >ip 

< P T >1tt 



(19) 



As < p\ >i p >< p\ >i7r, in the high density limit 
(Au-Au collisions k' > k) we expect a ratio larger than 1, 
as the experimental data show. 

As far as we approach the low density limit, the ratio 
decreases, becoming closer to 1 or even lower. 

A more detailed comparison with experimental data on 
Au-Au, d-Au collisions discussion of the forward rapidity 
region can be found in reference |24j . An overall reasonable 
agreement is obtained. It is very remarkable that such 
agreement is based on the universal behaviour of the pt 
distribution given by Equation (jl3|) . 



6 Bose- Einstein Correlations 

Most of the studies of two body Bose-Einstein (B.E) corre- 
lations have paid attention to the parameters R S ide, Rout, 
Rl and not to the strength of the correlation, defined by 
the chaoticity parameter 



C 2 (0,0) = A 



(20) 



Experimentally, due to Coulomb interference and to 
the necessary extrapolations there are many uncertainties 
in its evaluation, however, some trend of the dependence 
of A on the multiplicity can be established. First, SPS min- 
imum bias data for O-C, O-Cu, O-Ag and O-Au collision 
show that A decreases as the size of the target increases, 
from A = 0.92 up to A = 0.16. However for S-Pb and Pb- 
Pb central collisions, where the values of r\ are larger A is 
also larger, A ~ 0.5 — 0.7. 

This behaviour can easily explained in the percola- 
tion framework 1291. Each cluster can be considered as 




1.2 

n 



Fig. 5. A as a function of r\. Experimental points are for semi- 
central S-Pb collisions (filled triangle), 18% central Pb-Pb colli- 
sions (non filled circle) and 10% central Pb-Pb collisions (filled 
box) at SPS. 




Fig. 6. to as a function of r/. Experimental points are for S-Pb 
semicentral collisions and 9% central Pb-Pb collisions at SPS. 



a chaotic source A = 1, and the production of particles 
from several clusters can be seen as the superposition of 
chaotic sources. In this schemce, A = n s /n, being n s the 
number of pairs produced in the same cluster and Ut the 
total number of pairs. In this way, A is proportional to the 
inverse of the number of independent sources (clusters), 
therefore it decreases with the density up to the critical 
percolation value. ^From this critical value, it increases 
with density. This behaviour is shown in Fig. 5. 

Similar considerations can be done concerning the s- 
trength of the three body B-E correlations, u>. NA44 Col- 
laboration has obtained for S-Pb collisions, u> = 0.20 ± 
0.02 ± 0.19 and for central Pb-Pb collisions lj = 0.85 ± 
0.02 ± 0.21. STAR collaboration obtains for central Au- 
Au collisions values of u> close to 1. This sharp variation 
from uj = 0.2 to u> = 0.8 — 1 in a small range of rj is eas- 
ily explained in the framework of percolation of strings. 
Now u> becomes proportional to the inverse of the squared 
of the number of independent sources (clusters) what can 
acomodate stronger variation compared to the case of two 
body. In Fig. 6 our result is shown |80| . 
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7 Forward- Backward Correlations 

An useful observable to check the percolation approach is 
the forward-backward correlation measured by the quan- 
tity 
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Dfw =< tiftib > — < rip >< riB > 



(21) 



where ripim denotes the multiplicity in a forward (back- 
ward) rapidity interval. In order to eliminate the short 
range correlations, the forward and backward o intervals 
should be separated by at least one unit of rapidity. On 
general grounds, one can see that Dpw is proportional 
to the fluctuations on the number of independent sources, 
(or clusters in our case) At vei 7 l° w density, Dpw 

should be very small, increasing with the density up to a 
maximum related to the largest number of clusters. At 
very high density, there is essentially only are cluster and 
hence D FW becomes small again. 

There are some experimental data measuring the pa- 
rameter b, through 



< A*b >f- 



(22) 



where b = Dfb/DFF. The data on pp and pA show on 
increase of b with energy and density. Our prediction for 
high density is that b will decrease. Measurements of Dpw 
or 6 as a function of centrality would be welcome. 



8 Conclusion 

The percolation of partons and strings can describe rightly 
several observables, namely J/ip suppression, multiplic- 
ities, transverse momentum fluctuations, transverse mo- 
mentum distributions and B-E correlations. The behaviour 
of all of them has a common physical basis: the cluster- 
ing of color sources and the dependence of the number 
of cluster on the density. In this way, the threshold of 
J/tjj suppression, the maximum of transverse momentum 
fluctuations, the suppression of the Cronin effect and the 
turnover of the dependence of the strength of two and 
three body correlation with the energy are related to each 
other and all of them point out a percolation phase tran- 
sition. Another test of this transition is the measurements 
of forward-backward correlations and also the multiplicity 
distributions not discussed here |33| . 

Many of the results obtained in the framework of per- 
colation of strings are very similar to the one obtained in 
the color glass condensate (CGC). In particular, very sim- 
ilar scaling lows are obtained for the product and the ratio 
of the multiplicities and transverse momentum. For this 
reason, it is very tempting to identify the momentum Q s 
which established the scale in CGC with the correspond- 
ing arc in percolation of string. In this way 



_ k < p\ >i 



(23) 



The consequences of Eq. Q23JI are under study. 
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